Perfect fluids from high power sigma-models 



Radu Slobodeanu 



Department of Theoretical Physics and Mathematics, University of Bucharest, 
P.O. Box Mg-11, RO-077125 Bucharest- Magurele, Romania. 
E-mail: radualexandru.slobodeanu@g. unibuc.ro 

Abstract 

Certain solutions of a sextic sigma- model Lagrangian reminiscent of Skyrme model corre- 
' spond to perfect fluids with stiff matter equation of state. We analyse from a differential 

geometric perspective this correspondence extended to general barotropic fluids. 
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1 Introduction 



Relativistic hydrodynamics is an important theoretical tool in heavy-ion physics, astro- 
physics, and cosmology. In a first stage, it provides us with a good approximation in terms 
of perfect fluids for many physical systems, from quark- gluon plasma to gaseous interior of 
a collapsing star. A perfect fluid is defined [43J as having at each point a velocity v, such 
1^ ■ that an observer moving with this velocity sees the fluid around him as isotropic. For a 

. standard introduction in this subject see e.g. [TT t 

The fact that perfect fluid dynamics admits a Lagrangian formulation is widely known. 
In the approach taken in [71 [T^ [THt [TUl [5Dj , the Lagrangian is given in terms of a " dual" 
'sj" ' variable ip which is a submersion / field defined on the spacetime such that the fluid velocity 

four- vector spans its vertical foliation. In this (locally) dual picture the Euler equations for a 
I stiff fluid come by variational principle from a sextic Lagrangian in the derivatives of (p, and 

the conservation of energy along flow lines is automatically satisfied by construction. It is 
interesting to notice that this sextic Lagrangian is exactly the extension of the Skyrme model 
I discussed in [H [18]. In fact the equations for all perfect fluids with cosmological equation 

' of state come in a similar way by raising to a proper power this sextic Lagrangian. In this 

paper, we detail these constructions in the standard (coordinate-free) differential geometric 
setup and then we address the following topics: the geometric interpretation of the Euler 
equations for the fluid (Euler-Lagrange equations for the field), the effect of (bi)conformal 
change of metrics, the irrotational case, the shear-free case in relation with the harmonic 
morphisms theory and the coupling with gravity. We illustrate these issues with detailed 
examples that cover most of the presently known exact solutions for perfect fluid equations. 
The solutions are re-derived in this new context via the reduction by symmetry method 
that simplify Euler equations to an ODE. Their natural derivation and the fact that they 
automatically satisfy the conservation of energy along the flow are two primary advantages 
of the present approach that may further contribute to uncover other exact solutions yielding 
a new insight into the fluid dynamics |3.,2J. 

Our initial motivation was to bring the relativistic fluid problem in the intensively studied 
area of geometric variational methods and analysis of critical mappings between manifolds 
(mainly developed in the Riemannian case, see e.g. O US]), but the paper can also be 
connected with the kinetically driven inflation models [21 [31 [T2] or gravitating skyrmions 
[38|. 
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2 Perfect fluids 



Throughout the paper {M"^,g) will denote an m-dimensional Lorentz manifold. In par- 
ticular, the m-dimensional Minkowski space M"* = M^" is defined to be endowed 
with the metric of signature (l,m — 1) given in standard coordinates {xi,X2, ■■■,Xjn) by 
g = — dx^ + dxf + ... + dx'^_i (so Xm denote the time and Xi, i = 1, ...,m — 1 the spatial 
directions). By (m-dimensional) spacetime we shall mean a connected time-oriented Lorentz 
manifold. 

Recall that a smooth distribution V of dimension g on M is a smooth rank q subbundle 
of the tangent bundle TM. Assuming that V is non-degenerate, the metric g gives rise to a 
splitting TM = V ® and one can define the mean curvature of V by 



where V is the Levi-Civita connection, {cj} is a (local moving) frame for V, £i = g{ei, Ci) = 
±1 and denotes the orthogonal projection on "H of a vector field X on M. If V is closed 
under the Lie bracket, then it is called integrable distribution and, by Frobenius' theorem, 
it gives rise to a smooth foliation of M. The foliation / distribution V is said to be extremal 
if = (the terminology is inspired by the (local) "area" maximization property of zero 
mean curvature spacelike hypersurfaces, among all nearby hypersurfaces having the same 
boundary |19) . In the Riemannian case the same condition assures the minimality in a 
similar variational problem, see e.g. |42|). 

Now let us recall a slightly simplified picture of perfect fiuids, following |34[ Chap. 12]. 
Definition 2.1. A perfect fluid on a spacetime M is a triple {U,p,p) such that: 

(i) [/ is a timelike future-pointing unit vector field on M, called the flow vector field. 

(ii) /9,p : M — )• R are the mass (energy) density and pressure, respectively. 

(iii) divT = 0, where T is the stress-energy tensor given by: 



with Lo = U the 1-form metrically equivalent to U. 
Proposition 2.1. // (U,p,p) is a perfect fluid, then the following relations hold: 

(i) {p + p)divU + U{p) = conservation of energy along flow lines^ 

{ii) {p + p)'VuU + grad^p = (^Euler equations^, 
where grad^p (the spatial pressure gradient^ denotes the component orthogonal to U. 

In what follows we restrict ourself to barotropic fluids, i.e. we assume the equation of 
state (EoS): p = p{p). In this case, a third equation governing the fluid (the particle number 
conservation) decouples and therefore it is enough to work with the above deflnition. 



For reasons that will become clear below, we shall denote the foliation defined by U with 
V calling it vertical and its complementary with %, calling it horizontal distribution. Their 
mean curvatures are in this case given by: 




(2.1) 



1 = pg + {p + p)oj (g) oj 



(2.2) 



In this case it is useful to introduce the index of the fluid [llj by: f{p) := exp 





divU 



U. 



(2.3) 



m — 1 
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The fluid's equations above rewrite as foflows: 

^« = -grad^(ln(/')^); /^^ = grad^(ln /). (2.4) 

Example 2.1 (Radiation case). For instance, considering radiation's EoS, p = 3p and 
m = 4, we can resume fluid's equations in 

p^ + = grad(ln/), 

with / = p^/'*. The fact that + /^^ is of gradient type means, geometrically, that both V 
and % are extremal (i.e. = and = 0) with respect to the metric 'g = f'^g conformal 
to the spacetime metric. In particular, the flow lines are geodesies of cf. [11, Theorem 
10.1]. 



3 Perfect fluids from sigma-models 

3.1 The Cauchy-Green tensor 

Let if : {M"^,g) — )■ {N^,h) a map deflned on a spacetime with values in a Riemannian 
manifold. The pullback metric ip*h is a 2-covariant tensor field on M. If we define dip^ : 
TN — TM the formal adjoint of dip with respect to g and h in the usual way, then we can 
think of ip*h as an endomorphism = d<y9* o dip : TM — )• TM called the Cauchy-Green 
(strain) tensor of (p, analogously to the case of deformations in non-linear elasticity. It is 
self-adjoint by construction, but due to the Lorentzian signature of the metric, it may have 
complex eigenvalues. We shall eliminate this possibility by asking our maps to have timelike 
(so non-degenerate) vertical distribution V = kerdc^. 

For the convenience of the reader we present the following basic linear algebraic fact by 
specifying the dimensions to the case of interest for us (for general facts about self-adjoint 
linear operators on Lorentz vector spaces, see [Ml P- 243, 261]). 

Lemma 3.1. Let F : V ^ W be a linear operator of rank 3 from a Lorentz vector space 
V ^ R.f to an Euclidean vector space W ^ M^^. On V, consider the self-adjoint operator 
C := i?t o F, where F^ : W ^ V is defined by: {Fv,w) Euclid = {v, F^w) Minkowski, Vf G 
V,w£W. 

If {v,v) Minkowski < 0,Vf E ker i^, then the eigenvalues of ^ are (real and) non-negative 
and <t can be diagonalized in an orthonormal basis (of eigenvectors). 

Proof. Let J- = (/ij)i=i'2 3'^ be the associated (real) matrix of F, with respect to some 
orthonormal frames. Then the associated matrix of the adjoint -F* is the pseudo-transposed 
—fu — /21 — /31 \ 



matrix J^'^ 



fl2 /22 /32 

/i3 /23 /sa 

V /l4 /24 /34 / 



The associated matrix of the composition is a 



pseudo- symmetric matrix of the form J^* • = ^ ^ ^ J ' ^^^^^ a G M, ^ is a sym- 
metric matrix and zu a row 3- vector. Now let us apply a standard argument similar to the 
Euclidean case. Take Z = {zi)i<i<4 G a solution of the system S : (J^* • T — \I4)Z = 
and multiply it on the left with the diagonal matrix diag(— zi, Z2;^3;^4)- Summing the new 
equations provides us with a linear equation in A: A(— zi^i + Yli=2 ^i'^i)^ (real terms) = 0. 
So A will be real, unless we have some complex lightlike solution of S. But if it exists such 
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solution Z = X + iY, then {X, X) + (1", Y) = Q and X and Y are in the 1-dimensional kernel 
of F. So X and Y must be colinear and lightlike, that contradicts our hypothesis. 

The proof of the second statement is similar to the Euclidean case. ■ 

So, at a point x of M, if ker dipx is timelike, we saw that there exists an orthonormal 
basis of eigenvectors of for T^M. Extending [37^ Lemma 2.3] in this context, we get a 
local orthonormal frame of eigenvector fields, around any point where the timelike condition 
above is satisfied (i.e. (t^p can be consistently diagonalized). 

In general, let Ai, A2, A.^ and A^+i = ... = A^ = be the eigenvalues of £<^, 
where r := rank(d(/3) everywhere. Whenever they are all real, the corresponding elementary 
symmetric functions, ak{(p), will play a special role, due to the fact that 

ai{(p) = |d(/9p; aki^p) = \ dyjp; ... ; am{v) = det(y?*/i), 

where A^d(/5 denotes the induced map on fc-vectors and in the right-hand side of every 
identity we have the usual Hilbert Schmidt square norm induced by g and h (that may be 
negative, despite the squared notation). 



3.2 High power Lagrangians and their associated stress-energy tensor 

For a map if : (M, g) — t- {N, h) between (semi-) Riemannian manifolds let us recall 
some Lagrangians and their action integral £ (that will be called energy to recover the 
Riemannian geometric terminology). The associated stress-energy tensor S{ip) is defined by 
the first variation formula with respect to the domain metric variations: 

d_ 

du 



£{'P,9u) = / (5'((/?), 5g)i'g. 
«=o Jm 



(cf. |4J) A;-energy of ip: £k{^) = JMek{(p)i^g, where ek{(p) = l\dip\'', k>2. 
fc-stress-energy tensor: 

Sk{^):=eki^)g-\dip\''-''ip*h. (3.1) 

Adapting the proof in the Riemannian case, we easily get the identity: 

h{Tki^p),dip) = -dWSk{(p), (3.2) 

where Tfc((^) = trace^Vddc^l^^^d^j) is the k-tension field of 99. For submersions, (|3.2p 
tells us that Ski^p) is divergence free if and only if ip is k-harmonic, i.e. a stationary 
point of Ski') with respect to smooth variations of (p on any compact domain K 
{-j^\s=o£k{^s) = 0). A 2-harmonic map is called simply harmonic map or wave map 
if M is Lorentzian. 

(cf. [45J) cjfc-energy of v?: SakiV') = k hl'^k{^)^g, where ak{'p) = \ a'^ dc^p. A; > 1. 
(Tfc-stress-energy tensor: 

1 
2 



Scjki'P) = T,^k{'p)g - ^p*hoxk-i{^), (3.3) 



where Xq denotes the g*'* Newton tensor. Again, for submersions we can prove that 
Scrf.{(p) is divergence free if and only if ip is ak-critical, i.e. ■^\s=o£akiv^s) = on any 
compact domain K, or equivalently Ta^.{(p) = 0, where Tfj^[(p) = tracegV(d93o;i^;^_]^((^)) 
is called the Uk-tension field of <p>. 
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Notice that gO'il'/^) = ^2{y^) is the standard kinetic term and cr2{y^) coincides with the 
(self-interacting) Skyrme term (for fields with target N) . The stress-energy tensors of these 
two cases appear explicitly in [21] where one shows that they satisfy the dominant energy 
condition 0. 

It is worth to emphasize that the identities ()3.2p (consequences of the diffeomorphism 
invariance of the action) implies that the stress-energy tensor is always conserved along 
vertical directions (i.e. [divS'((/3)](V^) = 0,\/V £ Keidip) for any map (p, not necessarily 
for critical maps. In particular, for aa-energy we have the following general result, whose 
independent proof can be find in |27'. Prop. 3]. 

Lemma 3.2. Consider a submersion ip : {M^,g) — )• {N^,h) between (semi) Riemannian 
manifolds such that (T'i{ip) = det((/9*/i) >{) on an open set. Let U be an unit vertical vector 
field on this open set. Then the following identity holds: 

divU + [/(In ^/a^) = 0. (3.4) 



3.3 The (Ta-critical field / stiff" perfect ffuid correspondence 

In the following we identify a first case when the stress-energy tensor of a cr-model is in the 
perfect fluid form (see also [71 [12l [13] ) . This happens also for the stress-energy of a scalar 
field (kinetic or higher power) Lagrangian with potential, but in this case U is constrained 
to be irrotational, see e.g. [2| [15]. 

Proposition 3.1. Let (p : {M'^,g) — t- (iV^,/i) be a submersion from a Lorentzian manifold 
onto a Riemannian manifold. Suppose that the fibers are timelike and take U the unit 
vector field that spans the vertical space V = kerd(^. Denote by H = V"*" its complementary 
distribution and let g = g'^—g^ be the induced splitting of the metric tensor, with g^ = uj0uj, 
u}{X) := g{U,X), \/X G r(TM). Then the following statements hold good: 

(i) The eigenvalues of the Cauchy- Green tensor (tip with respect to g are all positive: 
Ai = \\,i = 1,2,3, A4 = 0. 

[ii) a ■i- stress- energy tensor of ip has the perfect fluid form 

-2S„,{^) = T=pg'^ + pg^, (3.5) 



with the stiff matter EoS: p = p = }^\\\\. 
{Hi) (p is CT^- critical if and only if 

-grad'^(ln A1A2A3) + = 0. (3.6) 

that is if and only if ([/, p, p) satisfies the relativistic Euler equation for the EoS 

p = p = xlxlxl. 

(iv) The following identity is true 

3^^ = grad^(lnAiA2A3). (3.7) 

and it is equivalent to the conservation of energy along flow lines for {U,p,p) with 
p = p = AfA^A^. 

(v) The fluid equations taken together are equivalent to: 

+ is of gradient type. 



^In our setting the stress-energy has the opposite sign with respect to the physics literature where it is 
defined as a variational derivative for the Lagrangian density relative to the inverse metric. 
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In particular after a biconformal rescaling of the metric both V and H become extremal. 

{vi) The equation ^3.b]) is invariant under biconformal changes of metric of the following 
type: 

5 = ?-V-?-V. (3.8) 

In particular, if Lp is a^-critical, then {c'U^cf' p,<f'p) is a stiff perfect fluid on {M,g). 

Proof, (i) Consequence of Lemma |3. 11 

{ii) Recall [35] that the second Newton tensor is given by: 

X2{^) = (T2{^p)IdTM -^ipo xiiy^y, Xi{^) = 2e2{ip)IdTM - 

where (J,^ = dip* o dip is the Cauchy-Green tensor. Let % = Hi ® 'H2 ® I-Lz be the splitting of 
the horizontal distribution in eigenspaces of C^. Take Xi G T{T-Li) and let us compute the 
(T3-stress-energy (j3.3p on it: 



s,, ((p) (Xi , y ) = , y ) - ^* h{x2 (v^) (Xi ) , y ) 

= \xlXl\lg{X^,Y) - ^*h{a2{ip)Xi - 2e2(95)C^(Xi) + <tl{Xi),Y) 

= IxlXlXlgiXuY) - ip*h{[XlXl + X^Xl + XlXj - (A? + + Xl)Xl + Xf\ X^Y) 

= XlXl (^^g^{X^,Y) - ip*h{XuY)^ = -^XlXlXlg^iXuY). 

A similar identity holds for the other two eigenspaces. 
Since for any V £ r(V) we have 

S^,iip){V,Y) = ^a;{ip)g{V,Y) - ^*h{x2mV),Y) = -iA?AiAi/(y, y), 

the conclusion follows. 

{Hi) Rewrite Euler equations from Proposition 2.1 in this context or simply take the 
divergence in (|3.5p . 

(iv) Apply Lemma [3. 2 1 Then interpret the mean curvature of H as we have done in ()2.3p . 
(f ) Under a general biconformal change of metric g = g^ — g"^ , the mean curvatures 
of horizontal and vertical distributions become: 

Jl-^ = ^[^,^ + gr^d^ (In g)]. ^ ^ ^ 

As in our case both fi^ and /x^ are of gradient type, they can be set to zero by such a 
biconformal change of metric. 

(vi) Using the second equation in (|3.9p . we can rewrite Equation (|3.6p with respect to a 
biconformally related metric g: 



- grad^(ln A1A2A3) + giad'^ (in c^'^ g) = 0. 



and the conclusion follows. 
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Remark 3.1. (a) The above result can be easily generalized for any dimension m of M, 
considering o"m-i-critical submersion with timelike one-dimensional fibers. In particular, 
the strongly coupled Faddeev-Niemi model [T7] on M'^ will have solutions interpretable as 
stiff perfect fluids in (2+1) dimensions. 

Notice that the condition for flbers to be timelike assures the dominant energy condition for 
the high power sigma-model considered (recall our different sign convention for the stress- 
energy tensor). For a more general statement see [44j . 

(b) A computation similar to (m) in the above proof, applied to the other possible energies 
eki^p) and ak{^p), shows us that the stress-energy tensor does not have the perfect fluid form 
(see also [2TJ). 

(c) Conversely, given a timelike unit vector fleld U, by locally integrating it becomes 
tangent to the fibres of a submersion onto some Riemannian manifold and the above result 
applies. So locally we are able to translate the problem for the triple {U,p,p) in terms of 
((/3, (A^, h)), once we have fixed a metric g on the spacetime M. 

(d) If A| = = A| = A^, then ip is called horizontally conformal [S] and the vertical 
unit vector U is shear-free. In this case (U,p,p) is a stiff perfect fluid if and only if (p 
is 6-harmonic (morphism). This condition translates into the following analogue of the 
fundamental equations for harmonic morphisms: 

{n — 6)grad^(ln A) + {m — n)fi^ = 0, 

where m = 4, n = 3 in our case (see the subsection 3.5 below). 

The particular case of semi-Riemannian submersions (i.e. A = 1) with extremal flbers 
{jjy = 0) will describe then fluids with constant pressure and mass density. 

(e) (Conformal change of the codomain metric). Let h = z7^/i, where z7 : — )• (0, oo) 
is a smooth function. Set v = 17 o (p. Then the eigenvalues of ^p*h are A? = (z^Aj)^. The 
equations ()3.7p leave unchanged, while ()3.6p become: 

-3grad^(lni/) - grad'^(ln A1A2A3) + /i^ = 0. 

Now let us illustrate the above discussion with two extensions of the well-known Hwa- 
Bjorken solution in (1+1) dimensions [9] (for a more elaborate treatment of the reduction 
techniques involved below, see [Fi Chap. 13]). 

Example 3.1 (Extended Hwa-Bjorken flow (I)). Let be the Euclidean space parametrized 
by spherical coordinates and regard as R x M^. Consider tp : — t- M^, an 5*0(3)- 
equivariant submersion of the form: 

(^X4, r(coss, sins-e*^)^ 1— )• a(x4, r) ^cos s, sins-e*^^, (3.10) 

depending on a function a : R x [0, 00) — [0, 00). Then the eigenvalues of £^ are: 

with the associated eigenvector flelds forming an orthonormal frame: 

^ 1 1^ da d da d \ ^ ^ X ^ U ^ ( da d da d\ 

^ Ai \dx/idx/i dr dr ) ^ ^ r ds^ ^ r sins 90' Ai \dr dx^ dx^dr J 

Notice that the condition for (f^)^ — to be positive on an open set translates 

both the condition of timelike flbers and the positivity of the flrst eigenvalue. 
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Writing U = cosh 0,{x4, r) ■ +sinh ^{x/^, r)--^, the equations p.6p reduce to one single 
PDE: 

\ox4, or J \ 0x4 or J 

Let us now restrict f to the interior of the forward hghtcone U = {x £ M^||x4| > |r|} 
and consider the ansatz 0(2:4, r) = /(C)j C •= r/x^. The common unit vertical vector for all 
maps in this ansatz is the (normalized) scale transformation: 

ir= ^4 _d_ r d_ 

./^l^dx4 ^^l^dr ^ ■ ^ 

Moreover, (fSTTTT) reduces to an ODE: [f'f + 2(/')V) C(l - O - 2/7^(1 + C^) = with 
the solution (C is an integration constant): 



The corresponding field is therefore crs-critical. As in this case A2A3 = {x\ — r^)~'^, 
the associated stiff perfect fluid is (C/,p = {x\ — r'^)~'^,p = p), cf. also [Hj. Recall that 
the energy conservation along flow lines is satisfled by construction as (j3.4p is true for any 
submersion. 

Since C/ is a standard conformal vector field, it is possible to construct a horizontally 
conformal submersion with fibers tangent to U. A way to do this is to firstly impose 
horizontal conformality (i.e. Af = A2) and then to employ Remark 3.1 (e) to render if as- 
critical. We illustrate this strategy for a co-rotational submersion to the 3-sphere ^p : M'^ — )• 
(S^, can) given by 

^X4, r(coss, sins-e*^)^ 1— )• ^cos /3(x4, r), sin/?(x4, r) (cos s, sins-e*' 

depending on a function /3 : R x (0,oo) — t- [0, vr]. Restricting again to an ansatz /3(x4,r) = 
/(r/x4) in the interior of the forward lightcone and asking (p to be horizontally conformal 
(i.e. Af = A2), we are leaded to an ODE with the explicit solution / = arcsin. In this case: 
Af = A2 = A| = X4 ^ and after a conformal rescale of the metric h = (cos~^ /3) • can, 
we shall obtain (via Remark 3.1 (e)) a as-critical map into a hemisphere, with the same 
associated perfect fluid. 

Let us comment some particular features of the above solution: 

• the flow is not accelerating. That is, the associated submersion have minimal fibers 
(VuU = 0) and the transversal "volume" Af A2A3 is constant in horizontal directions; 

• the flow vector U is shear-free, i.e. Cug^ — |div?7 ■ = 0; 

• the bulk viscosity of the fluid is divergence free, i.e. div(divf/ • (7^) = 0; 

• the heat-flow vector Q = giad^T +TV uU , for the standard temperature profile T = (x| — 
r^)~ 2 is zero. Together with the above two properties this means that this extension of Hwa- 
Bjorken flow is also a solution of the relativistic Navier-Stokes equations (see Appendix). 

Example 3.2 (Extended Hwa-Bjorken Row (II)). Let be the Euclidean space 
parametrized by cylindrical coordinates {x'^,x± ■ e^*^) and regard M"^ as M x M'^. Consider 
if : — )• R^, an S'0(2)-equivariant submersion of the following form: 



. X4, xs, x± ■ e 



^ (/3(x4,xs), xx-e*-^), (3.13) 
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depending on a function a : M x [0, oo) — )• M. Restricting ip to the future wedge X4 > |x3| and 
choosing the ansatz /3(a;3,X4) = /{x^/x^) we can proceed as in the previous example to find 
the solution /{x^/x^) = C arctanh(x3/x4) (C is an integration constant) for the reduced 
equation, solution that provide us with a era-critical submersion. The corresponding stiff 
perfect fluid is then given by: 

TT ^4 d X3 d 1 

U = ^^ =- h p = p=— 2 (3-14) 

^Jxl - xl OXi - x\ oxz x\ - x% 

Notice that the flow ()3.14p is no more shear-free and this solution does not trivially 
satisfy the relativistic Navier-Stokes equations as the previous one did. Nevertheless in [22] 
a correction of the energy density has been given such that the same flow proflle, with 
radiation EoS p = 3p, becomes an exact solution for the relativistic Navier-Stokes equations 
with X = C = ^iid non-zero shear viscosity coefficient. 

3.4 The crf-critical field / cosmological perfect fiuid correspondence 

In what follows we present the generalization of the above approach to the (linear) barotropic 
EoS: p = (7 — l)p, where 7 is a constant. When 1 < 7 < 2, this EoS is called cosmological 

m- 

Consider the following energy integral of a submersive field ip : {M'^,g) — )■ {N^,h) from 
a spacetime to some 3-dimensional Riemannian manifold: 



£^k{(f)= \A''d^rug, k>0. (3.15) 
Jm 



Theorem 3.1. The stress-energy tensor associated to i3.15\) is given by: 

S^.i^) = -^-^Xfxfxfg-^ - lx?Xfxfg\ (3.16) 
and the Euler- Lagrange equations associated to 13.15\) are equivalent to: 

-{2k - l)grad^(ln A1A2A3) + /i^ = 0. (3.17) 
Proof. Use again ()3.3p to compute 



du 



=0 "3 



'^a$(v2,9n)= / (T^^^U'P^a) ■ 9 - ka^ '^{^,9) ■ V*ho X2{v,g), Sg)iyg. 



M 



In order to conclude, perform a similar computation as in the proof of Proposition 13.11 (ii), 
then take the divergence of S^k{ip) on a horizontal vector. ■ 

Corohary 3.1. Let ip : {M'^,g) {N^,h) be a submersion with timelike fibers tangent 
to the the unit vector U. Then the variational equations { 3.17\ ) together with the identity 



represent, respectively, the relativistic Euler equation and the conservation of energy 
along flow lines for the perfect fluid {U,p,p) with p = Af'^Al^Al*^ and the cosmological EoS 
p={2k- l)p. 

In particular, we obtain a radiation EoS when k = 2/3 (when the corresponding - 
variational problem is conformally invariant), the cosmological constant EoS p = —p when 
A; = (identically constant Lagrangian) and the dust EoS p = when k = 1/2. 
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To illustrate the above results, simply notice that the submersions constructed in Ex- 
amples 13.11 and 13.21 are o"|-critical for all values of k. Accordingly, the perfect fluids cor- 
responding to Example 13.11 are {U,p = {x\ — r'^)~^^ , p = {2k — \)~^p) with U given by 

Let us rewrite these solutions in the customary coordinates for heavy-ion collisions: 



3^4~2;|, 7/ = arctanh— , x _\_ = \J x\ + x^, i;^> = arctan — , 

where the Minkowski metric reads: g = — dr^ + r'^drf + dx\ + x\d<f)'^. 

Example 3.3. The radiation type perfect fluid analogous to Example 13.11 is 

T d x_\_ d 1 
U= I =— + - — , p = ^p=— . (3.18) 



Corresponding to Example 13.21 we have the radiation type perfect fluid: 

U = ^, p = 3p=4- (3.19) 

A generalization of the perfect fluid (j3.19p is given in \12\ I23j . It corresponds to the 
^(^(Sj-equivariant, (T|-critical (for all k) submersion ip : dS^ x M — §^ x M, 

(sinhf), cosh £>(cos s, sin s • e"'^), r/) i— )• ((cos s, sins-e"'^), ??), 

where dS^ is the de Sitter 3-dimensional space, is the two-sphere and both domain and 
codomain are endowed with the product metric. The unit vertical vector \s U = d/dg and 
the eigenvalues of are Af = 1, A2 = A3 = cosh~^ g ; thus ip induces in particular the 
radiation type fluid {U, p = 3p = cosh^^/'^ g) corresponding to (T3~^''^-variational problem. 

k=2/3 

As dS^ X M is conformally equivalent to the Minkowski space and cjg -energy is confor- 

mally invariant, the above submersion produces a cxg"^'' '^-critical map (so a radiation type 
perfect fluid) on the (future wedge in the) Minkowski space. 

Remark 3.2. More general EoS can be treated analogously, by considering an energy of 
the type J^jF{\ d^\'^)ug (see also [3 Ha O [l6] ) . 

3.5 r-Harmonic morphisms / shear-free barotropic fluids correspondence 

In the following we extend Remark 13.11 (d) by showing that the r-harmonic morphisms 
are related to shear-free perfect fluids with various EoS. Recall [26j that a (r-)harmonic 
morphism is a (r-)harmonic and horizontally conformal map. For examples and more details 
on harmonic morphisms see [5l El [20] . 

Proposition 3.2. Let ip : {M'^^g) — )• (A^^,/i) he a submersion with timelike fibers from a 
spacetime to some Riemannian manifold and U an unit vertical timelike vector field. If ip 
is horizontally conformal, i.e. Af = A2 = A|, then 

{i) U is a shear-free vector field, 

(ii) ip is -critical if only if it is a 6k-harmonic morphism. If this is the case, then 
{U,p = {2k — 1)X^^, p = A^^) is a shear-free perfect fluid with the EoS p = {2k — l)p. 
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Proof, (i) As f is horizontally conformal, we have X'^g^ = (p*h and therefore Cu {X?g^^ = 
that gives us the result. 

{ii) For a harmonic morphism on a m-dimensional spacetime, the fundamental equations 
hold [24|: (n - 2)grad^(lnA) + (m - n)^i^ = 0. They can be generalized for r-harmonic 
morphisms as follows: 

(n — r)grad^(lnA) + (m — n)fj^ = 0, 

where m = 4 and n = 3 in our case. The conclusion easily follows. ■ 

For instance, if (/? is a harmonic morphism with dilation A^, then {U,p = — |A^,p = A^) 
is a shear-free perfect fluid with the EoS p = —3p (quintessence). Analogously, if (/9 is a 
4-harmonic morphism, then {U,p = ^A^,p = A^) is a shear-free perfect fluid with the EoS 
p = 3p (radiation). 

We emphasize the fact that locally there is a one to one correspondence between horizon- 
tally conformal maps and shear-free (vertical) vector fields, cf. Prop. 2.5.11]. According 
to the previous Proposition this correspondence holds also between r-harmonic morphisms 
and shear-free linear-barotropic fluid flows (general barotropic flows can be treated analo- 
gously by employing the notion of F-harmonic morphism [35]). Now recall that harmonic 
morphisms with one dimensional fibers defined on Riemannian manifolds are subject to var- 
ious classification results \10\ [36] . Similar results for r-harmonic morphisms in Lorentzian 
signature (with the condition of timelike fibers) would locally classify the associated shear- 
free perfect fiuids. The main result in |41] is an example of such classification, for radiative 
shear- free fluids that are coupled with gravitation (see next section). 

Example 3.4 (Morawetz flow). For the submersions defined by ()3.10p let us consider the 
(Tg^^-Euler-Lagrange equations ()3.17p that reduce to one single PDF analogous to (|3.1ip : 



dx4 3 dr J \dr 3 8x4 



V.2 



This time we impose the ansatz a(x4,r) = /(ln(^^: — )) defined on the forward lightcone 
interior. The vertical timelike vector of these submersions will be the (normalized) Morawetz 
vector field 



c = 4±4 « +Jfi!:^f (3.21) 

X4 — 0x4 — r^ or 

In this case Xf = (f)"^ /r"^ and A2 = A3 = /^/r^ and (|3.20p reduces again to an ODF: 
(In/'/^)' = —3. The corresponding solution a(x4,r) = 2 provides us with a hori- 

2 /3 

zontally conformal cjg -critical map, that is with a 4-harmonic morphism. Computing p 
= (AiA2A3)^/^ = A^ = (x\—r'^)~^ , we get the associated radiative perfect fluid (C/, p,p = ^p). 
This shear-free and accelerated (i.e. VuU 7^ 0) flow was pointed out in [33] where other 
exact solutions with VL(x4,r) = k ■ arctanh(r/x4), k G M can be found. Recall that the 
Morawetz vector field belongs to the family of special conformal transformations so this 
example may admit further generalizations. 

Example 3.5 (Skew projection's flow). Consider the cylinders 

Cf(q) = {x = (X4, (X3, xx-e'-^)) G I 4 < Q'^}, C\q) = {x = (xg, x±-e'^) G | < 
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in Minkowski and 3-Euclidean space, respectively. On the first consider the induced Lorentzian 
metric and on the second consider the deformed metric h = (1 — g^i^) • [dx| + dx^] +x]_d(tP'. 
Let 93 : Cf{q) — )■ C^{q) be the axially symmetric submersion given by: 

(x4, (X3, XX • e^*)) ^ (x3, XX • ei^'^-''"^)) . (3.22) 

We can easily verify that U = , ^ „ H — , ^ „ -St is the unit timelike vertical vector 

■y -y yji~q^x\ 9x4, ^i_g2^2^ 

field and the eigenvalues of the Cauchy-Green tensor of if are = A2 = A3 = 1 — g^x^ := A^ 
and A4 = 0. Now let us notice that the following identities hold: 

Vc/C/-grad^(lnA) = 0, divC/ + C/(ln A^) = (3.23) 

where the second is trivially satisfied as U is divergence free. Then (j3.23p tells us simulta- 
neously that [U, p = — p/3, p = A^) satifies the equations of a (shear-free) perfect fluid and 
that the associated submersion ip is a harmonic morphism. 

Example 3.6 (Radiation flow in Einstein universe). Consider the static spacetime 
(R X S^,g) parametrized as (i, coss • e"^^,sins • e"^^) so that the product Lorentzian metric 
reads g = -dt^ + ds^ + cos^ s d(/>| + sin^sd02- Take cji,cj2 < 1 and define the axially 
symmetric submersion (/^^^ : M x S'^ — )• by 

(t,coss • e''^\sins • e''^^) ^ (^oss • e'^^^-'^'^\ sin s ■ ei(<^2-<^2t))_ 

whose unit timelike vertical vector field is: 

1 f d d d 

U = — , = I — + Wl 77^ + UJ2 

I — ujf cos^ s — UJ2 sin^ s 



dt U<pi U(p2 

3 I, _ V-3 



Write the standard metric on S as /i = e* ® e* with e*'s given by: 
= ds, 

= (wf cos^ s + siv? s)^^/^(a;i cos^ s + uj2 sin^ s d<f>2) and 

= sin s cos s(a;f cos^ s + ^in^ ■s)~^/^(a;2 d(/)i — wi d02). 
Now consider the following non-standard metric on the codomain: 
h = A{s)~'^[A{s)e^ (g) + (g) + A{s)e^ e^], where A{s) := 1 - uf cos^ s - w| sin^ s. 
Then the eigenvalues for with respect to h are: Af = A^ = A3 = A{s)~^ := A^, A4 = 0. 

9/0 

We can check that 99^^ is a 4-harmonic morphism (so -critical submersion), and there- 
fore [/ is a shear- free perfect fluid with the energy density p = {l — ujf cos^ s — UJ2 sin^ 
and p = ^p. This solution has been analysed in [F, §3] where it is shown that it also satisfies 
the relativistic Navier-Stokes equations for an appropriate temperature profile. 



3.6 Irrotational flows 

In the following we characterize irrotational flows in terms of their associated field. 

Suppose that the submersive field with timelike fibers if : {M'^,g) — )• {N^,h) is a critical 
point for the functional £^k. Then, according to Corollary 3.1, the timelike unit vector U 
that spans the fibres is the fiow vector field of a perfect fluid. Let V = {l/f)U be the 
fundamental vector field of 93, where / is the index of the fluid. Denote its dual 1-form with 
T?. The next result and its proof are analogous to the Riemannian case, see [3 p. 341]. 

Proposition 3.3. The 2-form di? identically vanishes if and only if the horizontal distribu- 
tion is integrahle, i.e. U is orthogonal to a (spacelike) hypersurface. 
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So d-!? win be again called the integrability 2-form of ip. Notice that V and d-i? play a role 
similar to the dynamic A-velocity and the vorticity tensor |llj . 

An easy check shows us that the horizontal distribution is integrable for Examples 13.11 
and l3.2|, so U is irrotational in both cases. 



4 Coupling with gravity 

Let (M, g) be a 4-dimensional spacetime, with Ric*^ and Scal^ denoting its Ricci and scalar 
curvatures, respectively. We call perfect fluid coupled with gravity a triple {U, p, p) that solves 
the Einstein equation: 

Ric^^ - ^ScaF • g = aT, (4.1) 



where a is constant and T is the fluid stress-energy tensor given by (j2.2p . The existence of 
such solution impose severe restrictions on the curvature as shown by the following result 
(see e.g. [M]). 

Proposition 4.1. Let U be a future-pointing timelike unit vector field on the spacetime 
M . Then U is the flow vector field of a perfect fluid coupled with gravity if and only if 
Ric^(X, Y) is proportional to g{X, Y) and Ric*^(X, U) = for all X and Y orthogonal to 
U. In this case the energy density and the pressure of the fluid are: 



1 

P = - 

a 



Ric^^(;7,C/) + ^ScaP 



1 

^ 3a 



Ric^^ ([/,[/) - ^ScaF 



A standard case when the above conditions on Ricci tensor are satisfled is provided by 
the warped product (M = I x fN,g = —dt^ + • /i) where ( A^, h) is a (Riemannian) Einstein 
manifold and / a function deflned on the real interval /. 

Turning to the variational picture, let f : (M^^^g) — t- {N^,h) be a smooth map from a 
Lorentzian manifold to a Riemannian 3-manifold. Consider a common energy £{g,ip) for 
the field ip coupled to gravity given as 

£{9, = £grav{g) " aSfieldig, (4.2) 

where £grav{g) = Jm Scal*^fg, £fieidig,f) = Jm I d(^p*'Vg and a is a coupling constant. 
For a similar coupling of the standard kinetic term with gravity, see e.g. |25 [ I28 [ [M]. 

Then the stationary points of £{g,ip) with respect to smooth variations of g and 93, on 
any compact domain K C M, are characterized by: 

r Ric*^ - iScaF • 5 = aS.(^); 

\ ' (4.3) 

I r^li^) =0, 

where S^k{ip) and T^k{ip) are the stress-energy tensor and Euler-Lagrange operator associ- 
ated to £ field- 

According to Theorem 13.11 if a submersion with timelike fibers satisfies the first equa- 
tion in (j4.3p then it will provide us with a (linear) barotropic fiuid coupled to gravitation. 
Notice that the second equation in (j4.3p will be a consequence of the first one. 

Let us give a simple example. More examples of (73 -critical submersions coupled with 
gravity can be obtained from the perfect fiuid solutions of the Einstein equations discussed 
in HOl. 
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Example 4.1. Let us consider the Robertson- Walker spacetime M^j = I Xf N, where 
(A^, h) is a 3-dimensionaI, connected Riemannian manifold of constant curvature k. Then 
U = d/dt is the flow vector of a barotropic perfect fluid coupled with gravity, with p and 
p depending only on the time i G / |34] (in fact it solves the full relativistic Navier-Stokes 
equations as the dissipative terms of the stress-energy tensor vanish in this case). 

If we take f{t) = at + h an affine function, then p = ?>{a? + k){at + and the EoS for 
this fluid is: p = —^p. The associated o"3~^''^-critical submersion is the projection on the 
second factor tt2 : (M^j, -dt^ + fit)h) {N, 3(0^ + k)-h); it is, in particular, a harmonic 
morphism (horizontally homothetic with totally geodesic fibers) and a solution of ()4.3p for 
k = 1/3 and a = —2. 

If we take f{t) = \Jlat — kt'^, a > (defined only for finite time when A; > 0), then p = 
3a^/(t)~^ and the EoS for this fluid is: p = ^p. The associated (T3~^^^-critical submersion 
IS again the projection 7r2 : (M^j, -dt^ + P{t)h) {N, a^/3 ■ h); it is also a 4-harmonic 
morphism and a solution of ()4.3p for k = 2/3 and a = —2. 



5 Appendix: Thermodynamic degrees of freedom 

Recall that generally a relativistic fluid is characterized by five thermodynamic scalars that 
are subject to the first law of thermodynamics [Gihhs equation): 

dp = ^ dn + nTds. (5.1) 
n 

Here n is the particle number density with the associated conservation law d\v{nU) = 0, 

s is the specific entropy and T is the temperature. Gibbs equation shows that in general 

two of these five scalars are needed as independent variables. For example, taking n and 

s as independent, the remaining thermodynamical scalars are p{n,s), p{n,s), T{n,s), and 

given any one of these, say p = p{n, s), the others will be determined hy p = —p + and 
rp _ 1 dp 

n ds' 

In our context n = A1A2A3 (so that the conservation law div(nC/) = is automatically 
satisfied according to Lemma l3.2p and p = F{n^), so p = —F{n'^) + 2n^-F'(n^). Therefore 
a fluid that comes from a sigma model will always have a barotropic EoS for the energy 
density and the pressure. It can be easily shown that in this case the fluid is isoentropic 
(i.e. s is constant), while a general perfect fluid is only adiabatic (i.e. U{s) = 0). It would 
be interesting to find an extension of the action principle used in this paper that account 
for entropy production (at least by shear and bulk viscosities). 

Recall that the dissipative fluids (or imperfect fluids) dynamics is modelled by the rela- 
tivistic Navier-Stokes equations divT = 0, where T is the modified stress-energy tensor 
p. 567], [ISl p. 56]) given by: 



pg+ {p + p)uj (SlUJ -rj {^Cug^ - ^divf/ • - ^(6* fX" w + w (g) 6*) - (divU ■ 5^. 

where 6 = (grad^T + TVijU)^ is the dual of the heat-flow vector and rj, X) C are the 
shear viscosity, heat conduction and bulk viscosity coefficients, respectively (all are scalars 
that may depend on the temperature and the mean free time). The above stress-energy 
choice is motivated by the fact that together with (|5.ip it assures a positive rate of entropy 
generation: div{nsU — y^") > {the second law of thermodynamics). 
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